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1. Introduction. 



In this work we present a study of the Ermakov-Lewis invariants that are 
related to some linear differential equations of second order and one variable 
which are of much interest in many areas of physics. In particular we shall 
study in some detail the application of the Ermakov-Lewis formalism to 
several simple Hamiltonian models of "quantum" cosmology. There is also a 
formal application to the physical optics of waveguides. 

In 1880, Ermakov [1]] published excerpts of his course on mathematical 
analysis where he described a relationship between linear differential equa- 
tions of second order and a particular type of nonlinear equation. At the 
beginning of the thirties, Milne @ developed a method quite similar to the 
WKB technique where the same nonlinear equation found by Ermakov oc- 
curred, and applied it successfully to several model problems in quantum 
mechanics. Further, in 1950, the solution to this nonlinear differential equa- 
tion has been given by Pinney 0]. 

On the other hand, within the study of the adiabatic invariants at the 
end of the fifties, a number of powerful perturbative methods in the phase 
space have been developed . In particular, Kruskal ||| introduced a certain 
type of canonical variables which had the merit of considerably simplifying 
the mathematical approach and of clarifying some quasi-invariant structures 
of the phase space. Kruskal's results have been used by Lewis P to prove 
that the adiabatic invariant found by Kruskal is in fact a true invariant. 
Lewis applied it to the well-known problem of the harmonic oscillator of 
time-dependent frequency. Moreover, Lewis and Riesenfeld JJJ proceeded 
to quantize the invariant, although the physical interpretation was still not 
clear even at the classical level. In other words, a constant of motion without 
meaning was available. 

In a subsequent work of Eliezer and Gray 0, an elementary physical 
interpretation was achieved in terms of the angular momentum of an aux- 
iliary two-dimensional motion. Even though this interpretation is not fully 
satisfactory in the general case, it is the clearest at the moment. 

Presently, the Ermakov-Lewis dynamical invariants are more and more 
in use for many different time-dependent problems whose Hamiltonian is a 
quadratic form in the canonical coordinates. 
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2. The method of Ermakov. 

The Ukrainian mathematician V. Ermakov was the first to notice that some 
nonlinear differential equations are related in a simple and definite way with 
the second order linear differential equations. Ermakov gave as an example 
the so-called Ermakov system for which he formulated the following theorem. 

Theorem IE. If an integral of the equation 

is known, we can find an integral of the equation 

d 2 z a 

where a is some constant. 

Eliminating M from these equations one gets 



d ( dz dy\ ay 
dx \ dx dx ) z 3 



Multiplying both sides by 



the last equation turns into 



2 / <^±_ z ^ 
V dx dx , 



d ( dz dy\ 2 lay d fy 
dx \ dx dx J z dx \z , 

Multiplying now by dx and integrating both sides we get 



^dz z dyV = c ay 
dx dx I z 2 



,2 



If yi and yi are two particular solutions of the equation (1), substituting 
them by y in the latter equation we get two integrals of of the equation (2) 



— - z—Y = C - — * 

dx dx l z 
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dz dy 2 V 

V2-T- - Z— =C 2 - 

dx dx J 



Eliminating dz/dx from these equations we get a general first integral of 
(2). One should note that the Ermakov system coincide with the problem 
of the two-dimensional parametric oscillator (as we shall see in chapter 6). 
Moreover, the proof of the theorem gives an exact method to solve this 
important dynamical problem. 

The general first integral of equation (2) can be also obtained as follows. 
Getting dx from equation (3): 



This is the general first integral of equation (2), where C 3 is the constant 
of the last integration. For y is enough to take any particular integral of 
equation (1). 

As a corollary of the previous theorem we can say that 

Corollary lEc. If a particular solution of (2) is known, we can find the general 
solution of equation (1). 

Since it is sufficient to find particular solutions of (1), we can take C = 
in equation (3). Thus we get: 



ydz — zdy 



\JC — ay 2 / z 2 



Dividing both sides by y 2 we get the form 




Multiplying by C and integrating both sides we get: 





and therefore 



dy dz dx 

y z 
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Integrating both sides 

dx 



logy = \ogz± y^a ( 

J z 



which results in 

y = zexp(± v /z a ^ ^). 

Taking the plus sign first and the minus sign next we get two particular 
solutions of equation (1). 

A generalization of the theorem has been given by the same Ermakov. 

Theorem 2E. Ifp is some known function of x and f is any other arbitrary given 
function, then the general solution of the equation 

p d 2 y _ y (Pp_ = }_Jv_ 
dx 2 dx 2 p 2 \p 

can be found by quadratures. 

Multiplying the equation by 



\ \JjJU \JjJU ) 



one gets the following form 



d ( V ^L - = 2f (-\d (-} 

\ dx dx J \P ) \P ) 

Integrating both sides and defining for simplicity reasons 

2 J f(z)dz = <p(z), 

we get 

/ dy _ dp\ 2 = fy\ +c 
\ dx dx J \p J 

This is the expression for a first integral of the equation. Thus, for dx we 
have: 

pdy - ydp 

dx = 



Dividing by p 2 and integrating both sides we find: 

d* c f "(J) 

— + W - 



P 2 



This is the general integral of the equation. 

A particular case is when p = x. Then, the differential equation can be 
written as 

3 d2 y t (V 



dx 2 \x 

3. The method of Milne. 

In 1930, Milne introduced a method to solve the Schrodinger equation 
taking into account the basic oscillatory structure of the wave function. 
This method has been one of the first in the class of the so-called phase- 
amplitude procedures, which allow to get sufficiently exact solutions for the 
one-dimensional Schrodinger equation at any energy and are used to locate 
resonances. 

Let us consider the one-dimensional Schrodinger equation 

^ + * a (*W = (1) 

where k 2 (x) is the local wave vector 

k 2 {x) =2fi[E-V(x)] . (2) 

Milne proposed to write the wave function as a variable amplitude multiplied 
by the sinus of a variable phase, i.e., 

/2u\ 1/2 

^( x ) = ( -f ) P( x ) sin(0(x) + 7 ) (3) 

where \x is the mass parameter of the problem at hand, E is the total energy 
of the system, 7 is a constant phase, and V(x) is the potential energy. In 
the original method, (3 and are real and (3 > 0. Substituting the previous 
expression for if> in the wave equation and solving for d<p/dx one gets 



d 2 (3 



+ k 2 (x)(3=-, (4) 



dx 2 x J ^ (3 3 
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dx I3 2 ' w 

As one can see, the equation for Milne's amplitude coincides with the non- 
lineal equation found by Ermakov, now known as Pinney's equation. 



4. Pinney's result. 

In a brief note Pinney was the first to give without proof (claimed to be 
trivial) the connection between the solutions of the equation for the para- 
metric oscillator and the nonlinear equation (known as Pinney's equation or 
Pinney-Milne equation). 

y"(x)+p(x)y(x) + ^(x) = (1) 

for C = constant and p(x) given. The general solution for which y(x ) = y , 
y'M = V'o is 

yp (x) = \u 2 {x) - CW- 2 V 2 (x)] 1/2 , (2) 
where U and V are solutions of the linear equation 



y"{x) +p(x)y(x) = , (3) 

for which U(xo) = yo, U'(x ) = y' Q ; V(x ) = 0, V'(xq) ^ y' and W is the 
Wronskian W = UV — U'V = constant ^ and one takes the square root 
in (2) for that solution which at x has the value y 

The proof is very simple as follows. From eq. (2) we get y' P = y^}{UlJ — 
CW~ 2 VV) andy'p = -y P \UU -CW- 2 VV) 2 +y P 1 {U 2 -CW' 2 V 2 )-p{x)y P . 
From here, explicitly calculating y P + p(x)yp one gets —Cyp 3 and therefore 
Pinney's equation. 



5. Lewis' results. 

In 1967 Lewis considered parametric Hamiltonians of the standard form: 

H L = (l/2e)[p 2 + n 2 (t)q 2 ] , (1) 
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If is real, the motion of the classical system whoose Hamiltonian is given 
by eq. (|l|) is oscillatory with an arbitrary high frequency when e goes to zero. 
Corresponding to this, there are asymptotic series in positive powers of e, 
whose partial sums are the adiabatic invariants of the system; the leading 
term of the series is eH/Q. In the problem of the charged particle the adi- 
abatic invariant is the series of the magnetic moment. Lewis's results came 
out from a direct application of the asymptotic theory of Kruskal (1962) to 
the classical system described by Hl with real Q. Lewis found that Kruskal's 
theory could be applied in exact form. As a consequence, an exact invariant, 
which is precisely the Ermakov-Lewis invariant, has been found as a special 
case of Kruskal's adiabatic invariant. Although Q was originally supposed 
to be real, the final results hold for complex Q as well. Moreover, the exact 
invariant is a constant of motion of the quantum system whose Hamiltonian 
is given by the quantum version of eq. (|T]). 

The classical case. 

Let us take a real Q. In order to correctly apply Kruskal's theory, it is 
necessary to write the equations of motion as for an autonomous system of 
first order so that all solutions be periodic in the independent variables for 
e — > 0. This can be achieved by means of a new independent variable s 
defined as s = t/e and formally considering t as a dependent variable. The 
resulting system of equations is 

dq/ds = p , 
dp/ds = -tt 2 (t)q , 

dt/ds = e (2) 

Since t is now a dependent variable, this system is autonomous. In the 
limit e — > 0, the solution of the last equation is t — constant, and therefore 
the other two equations correspond to a harmonic oscillator of constant fre- 
quency. Since fl is real, the dependent variables are periodic in s of period 
2n/Q(t) in the limit e — 0, and the system of equations has the form required 
by Kruskal's asymptotic theory. A central characteristic of the latter theory 
is a transformation from the variables (q,p,t) to the so-called "nice vari- 
ables" (zi, 22 j V 9 )- The latter are chosen in such a way that a two-parameter 
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family of closed curves in the space (q, p, t) can be defined by the conditions 
Z\ = constant and z<i = constant. These closed curves are called rings. The 
variable ip is a variable angle which is defined in such a way as to change 
by 27r if any of the rings is covered once. The rings have the important fea- 
ture that all the family can be mapped to itself if on each ring s is changed 
according to eqs. (0). In the general theory, the transformation from the 
variables (q,p, t) to the variables (zi, z 2 , (p) is defined as an asymptotic series 
in positive powers of e, and a general prescription is given to determine the 
transformation order by order. As a matter of fact, Lewis has shown one 
possible explicit form for this transformation in terms of the variables q, p 
and Pinney's function p(t) . Moreover, the inverse transformation can also be 
obtained in explicit form. 

For the parametric oscillator problem, the rings are to be found in the 
t= constant planes. It is this property that allows the usage of the rings for 
defining the exact invariant I as the action integral 

I = <t> pdq . (3) 

J ring 

By doing explicitly the integral of I as an integral from to 2tt over the 
variable ip (see the Appendix), one gets 

I= l -[(q 2 /p 2 ) + (px-epq) 2 } , (4) 

where p satisfies Pinney's equation 

e 2 p + Q 2 (t)p - 1/p 3 = , (5) 

and the point denotes differentiation with respect to t. The function p can be 
taken as any particular solution of eq. Although Q was supposed to be 
real, / is an invariant even for complex Q. It is easy to check that dl/dt = 
for the general case of complex Q by performing the derivation of eq. @, 
using eqs. (0) to eliminate dq/dt and dp/dt, and eq. (^) to eliminate p. 

It might appear that the problem of solving the system of linear equations 
given by eqs. (§) has been merely replaced by the problem of solving the 
nonlinear eq. (|5]). This is however not so. First, any particular solution of 
eq. (|^) can be used in the formula of I with all the initial conditions for the 
eqs. (0). For the numerical work it is sufficient to find a particular solution 
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for p. Second, the exact invariant has a simple and explicit dependence on 
the dynamical variables p and q. Third, taking into account the fact that 
e 2 is a factor for p in eq. one can obtain directly a particular solution 
for p as a series of positive powers of e 2 . If O is real and the leading term 
of the series is taken as f2 _1//2 , then the corresponding series solution is just 
the adiabatic invariant expressed as a series. It is interesting to speculate if 
in practice it is more useful to calculate / by means of the solution written 
as a truncated series of p or by the corresponding expression in series for / 
truncated at the same power of e. Forth, one can also solve eq. (f|) to get p 
as a power series in 1/e 2 in terms of integrales. Finally, with the result of 
eqs. (|J) and (||), it is possible to get a better understanding of the nature of 
Kruskal's adiabatic invariant. Some progress in this regard can be found in 
the following general discussion on I and p. 

By adding a constant factor, the invariant / of eq. (^) is the most general 
quadratic invariant of the system whose Hamiltonian given by eq. is also 
a homogeneous quadratic form in p and q. This can be seen by writing the 
invariant in terms of two linear independent solutions, f(t) and g(t) of the 
parametric equation. If we write the generalized form of / 

I = 5 2 [p- 2 q 2 + (pp-epq) 2 } , (6) 

where 5 is an arbitrary constant, and compare this form with that in terms 
of f(t) y g{t), then we can infer that the two invariants are identical if p is 
given by 



2 r?2 -A 2 B 2 



A 2 „ B 



P = li(ea)- 1 —g 2 + —f 2 + 2 l2 —-(ew) 2 fg , (7) 



1/2 



1/2 



where A and B are arbitrary constants, while the constants a, 71 and 72 are 
defined by 

w = fg - gf , 71 = ±1 , 72 = ±1 . (8) 

Since there are two arbitrary constants, this formula for p gives the general 
solution of eq. (|5p expressed in terms of / and g. Using this formula we can 
build p explicitly for any Q for which the eqs. (0) can be solved in an exact 
manner. By constructing p in this way for special cases, we can infer that 
the expansion of p in a series of positive powers of e 2 is at least sometimes 
convergent. For example, if VI = bt~ 2n ^ 2n+1 \ where b is a constant and n is 
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any integer, the series expansion is a polynomial in e 2 , and consequently it 
is convergent with an infinite radius of convergence. 

Once we have the explicit form of Kruskal's invariant, it is possible to 
find a canonical transformation for which the new momentum is the invariant 
itself. If we denote the new coordinate by qi, the conjugated momentum by 
Pi, and the generating function by F, then the results can be written as 

q 1 = -tan~ 1 [pV<? - epp] , 

Pi = ^[p~V + (pP~ep<?) 2 ] , 
F = \tp' X pq 2 ± P~ X ^Pi - p- 2 q 2 ) 2 ± Pl sin" V V(2pi) 1/2 ] + (n + ^)n Pl 
~ < sin _1 [p _1 g/(2pi) 1/2 ] < | ,n = integer^) , 



dF OF 
qi = - 

dF 1 



P = ~R > 91 o 
oq opi 



H new = H+^ = ^p- 2 Pl . (9) 

In the expression for F the upper or lower signs are taken according to 
p — ep _1 pq is greater or less than 0. One can see that q\ is a cyclic variable 
in the new Hamiltonian, as it should be if p\ = I can be an exact invariant. 

Moreover, Lewis noticed that the second order differential equation for q, 
namely e 2 d 2 q/dt 2 + Q 2 (t)q = 0, is of the same form as the ID Schrodinger 
equation, if t is considered as the spatial coordinate and q is taken as the 
wave function. For bound states, Q is imaginary whereas for the continuous 
spectrum Q is real. Thus, the / invariant is a relationship between the wave 
function and its first derivative HI. 



The quantum case. 

Let us consider the quantum system with the same Hamiltonian H^, where 
q and p should fulfill now the commutation relations 

[q,p]=ih. (10) 
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We shall take p as real, which is possible if Q 2 is real. Using the commutation 
relationships and the equation for p it is easy to show that I is a quantum 
constant of motion, i.e., it can be an observble since it satisfies 

^ = #+4-[/,tf]=o. (ii) 

dt dt %h v J y ' 

It follows that / has eigenfunctions whose eigenvalues are time- dependent. 
The eigenfunctions and eigenvalues of I can be found by a method which 
is similar to that used by Dirac to find the eigenfunctions and eigenvalues 
of the harmonic oscillator Hamiltonian. First, we introduce the raising and 
lowering operators, a) and a, defined by 

= (l/y/2)\p- 1 q-i(pp-epq)\ , 

a = (l/V2)[p~ 1 q + t(pp-epq)] . (12) 
These operators fulfill the relationships 

[a, a'] = h , 

atf = I+^h. (13) 

The operator a acting on an eigenfunction of I gives rise to an eigenfunction 
of I whose eigenvalue is less by h with respect to the initial eigenvalue. 
Similarly, at acting on an eigenfunction of I raises the eigenvalue by H. Once 
these properties are settled, the normalization of the eigenfunctions of I can 
be used to prove that the eigenvalues of I are (n + ^)h, where n is or 
a positive integer. If |n) denotes the normalized eigenfunction of / whose 
eigenvalue is (n + j)h, we can express the relationhip between \n + 1) and 
\n) as follows 

\n + l) = [(n + l)h]-^ 2 ^\n) . (14) 

The condition that determines the eigenstate whose eigenvalue is is given 
by 

a|0)=0. (15) 

Using these results one can calculate the expectation value of the Hamiltonian 
in an eigenstate |n). The result is 

(n\H\n) = (l/2e)(p- 2 + Q 2 p 2 + t 2 p 2 )(n + \)h . (16) 
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It is interesting to note that the expectation values of H are equally spaced 
at each moment and that the lowest value is obtained for n = 0, i.e., we have 
an exact counterpart of the harmonic oscillator. As a matter of fact, we can 
obtain the harmonic oscillator results if Q is taken real and constant with 
p = fi -1 / 2 , which gives I = eH/VL. 



6. The interpretation of Eliezer and Gray. 



The harmonic linear motion corresponding to the ID parametric oscillator 
equation can be seen as the projection of a 2D motion of a particle driven 
by the same law of force. Thus, the 2D auxiliary motion is described by the 
equation 

^ + fi 2 (t)f=0 (1) 

where r is expressed in Cartesian coordinates (x, y). Using polar coordinates 
(p, 9) where p = \r\, x — pcos9, y = psin.9. The radial and transversal 
motions are described now by the equations 



p - p9 2 + n 2 P = o 

p dt 



(2) 
(3) 



Integrating eq. (H) 

p 2 9 = h (4) 

where h is the angular momentum, which is constant. Substituting in eq. (|^) 
one gets a Pinney equation of the form: 



2 h 2 

p + n 2 P = — 

p 3 



The invariant / corresponding to the eq. (§) is: 

.2 + (PP ~ X f>f 



1 = 1 

2 



h 2 x 2 



and with the substitutions x = p cos 9 and p = x one gets: 



1 



h 2 cos 2 9 + h 2 sin 2 9 



-h 2 



(5) 



(6) 



(7) 
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Thus, the constancy of / is equivalent to the constancy of the auxiliary 
angular momentum. 

In the elementary classical mechanics, the study of the simple ID har- 
monic oscillator is often made as the projection of the uniform circular motion 
on one of its diameters. The auxiliary motion introduced by Eliezer and Gray 
is just a generalization of this elementary procedure to more general laws of 
force. 

The connection between the solutions of the parametric oscillator linear 
equation and Pinney's solution is given by the following theorem. 

Theorem 1EG. If y\ and y 2 are linear independent solutions of the equation 

^ + Q(x)y = (8) 

and W is the Wronskian j/iy 2 — 2/22/1 (which, according to Abel's theorem is constant), then 
the general solution of 

d 2 y X 

where X is a constant, can be written as follows 

VP = (A Vl 2 + By 2 2 + 2Cy lV2 ) 1/2 (10) 
where A, B and C are constants such that 

AB -° 2 = W 2 



However, it is necessary that these constants be consistent with the initial 
conditions of the motion. If x\ (t) and x 2 (t) are linear independent paramet- 
ric solutions of initial conditions x\ (0) = 1, x\ (0) = 0, X2 (0) = 0, x 2 (0) = 1, 
the general parametric solution can be written as 

x (t) = axx (t) + f3x 2 (t) (12) 

where a and (3 are arbitrary constants that are related to the initial condi- 
tions of the motion by x (0) = a and x (0) = j3. The corresponding initial 
conditions for p and p are obtained from x = pcosO, x = pcosO — p9 sin 9, 
where 9 (0) = gives p (0) = a and p (0) = 0. Using ([TOD we get 

1 

(13) 



p(t)= (axi + (3x 2 ) 2 + y^2~ x 2 
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as the solution of (|5]) corresponding to the general parametric solution (|T2|). 

Moreover, we have 

p cos 9 = axi + (3x 2 (14) 

psin # = ^ (15) 

a 

The previous considerations can be extended to systems whose equations 
of motion are of the form 

§ + PW§ + QW^0. (16, 



The I invariant is now 

, h 2 x 2 



P 2 



+ {px - ppf exp (2 jf* P (t) dtj (17) 



where p is any solution of 

^ + P( t )| + (i)p= ^exp(- 2 {p( t )*) (18) 

The theorem that connects the solutions of (|16|) with those of ( |lq) (with a 
change of notation) can be formulated in the following way. 

Theorem 2EG. If y\ (x) and 2/2 (x) are two linear independent solutions of 

g + P(x)| + Q(x )y = (19) 



the general solution of 
ft 



^| + P(x)^ + Q(x) 2/ = A cxp f_ 2 / p( t )dt 
ax^ ax V \ J 



(20) 



can oe written down as 

y = [A Vl 2 + Py 2 2 + 2C7y iy2 ) 1/2 (21) 
where A and B are arbitrary constants, and 

A 



AB-C^ = ^exp(-2 / P(i)dt) (22) 
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7. The connection between the Ermakov in- 
variant and Nother's theorem. 



In 1978, Leach fllOfl found the Ermakov- Lewis invariant for the aforemen- 



tioned parametric equation with first derivative 

x + g(t)x + u 2 (t)x = , (1) 

by making use of a time- dependent canonical transformation leading to a 
constant new Hamiltonian. That transformation belonged to a symplec- 
tic group and has been put forth without details. In the same 1978 year, 
Lutzky [11 proved that the invariant could be obtained starting from a di- 



rect application of Noether's theorem (1918). This famous theorem makes 
a connection between the conserved quantities of a Lagrangian system with 
the group of symmetries that preserves the action as an invariant functional. 
Moreover, Lutzky discussed the relationships between the solutions of the 
parametric equation of motion and Pinney's solution and commented on the 
great potential of the method for solving the nonlinear equations. 

For the parametric equation without first derivative Lutzky used the fol- 
lowing formulation of Noether's theorem. 

Theorem NL. Let G be the one-parameter Lie group generated by 

such that the action functional J L(x,x,t)dt is left invariant under G. Then 

dL dL ,. . • dL 

where f = f(x,t), and 



BL dL ,, ■. dL • • 
^+»75T + (*-*0-5T+^=/. (2) 



• <9£ .<9£ . _dn dn ■ _ df_ . df_ 

at ox at ox at ox 



Moreover, a constant of motion of the system is given by 



dL 

* = {tix-n)— -ZL + f . (3) 
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The Lagrangian L = \{x 2 — oj 2 x 2 ) gives the equations of motion of the 
parametric oscillating type; substituting this Lagrangian in (H) and equating 
to zero the coefficients of the corresponding powers of x, on gets a set of 
equations for £, n, /. Next, it is easy to prove that they imply that £ is a 
function only of t and fulfills 

£ + Auo 2 i = . (4) 

The following results are easy to get 

n(x,t) = ^£x + ij(t) , 

f(x,t) = -for 3 + ipx + C , i> + u 2 tf) = . 

Choosing C = 0, ip = 0, and substituting these values in (||), one can find 
that 

*=\(& + [t<J 2 + ±£\x t -ixx) (5) 

is a conserved quantity for the parametric undamped oscillatory motion if £ 
satisfies ([|). Notice that (|j) has the first integral 

d - \i 2 + = c x . (6) 

If we choose £ = p 2 in ([5|) and (|]), with Ci = 1, we get that $ is the Ermakov- 
Lewis invariant. If the formula for the latter is considered as a differential 
equation for x, then it is easy to solve in the variable x/p; the result can be 
written in the form 

x = p[A cos (f> + B sin (ft] , <f> = <f>(t) , (7) 

where = 1/p 2 and A and B are arbitrary constants. Thus, the general 
parametric solution can be found if a particular solution of Pinney's equation 
is known. 

Consider now the Ermakov-Lewis invariant as a conserved quantity for 
Pinney's equation; this is possible if x fulfills the parametric equation of 
motion. This standpoint is interesting because it provides an example of how 
to use Noether's theorem to change a problem of solving nonlinear equations 
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into an equivalent problem of solving linear equations. Thus, if we take as 
our initial task to solve Pinney's equation, we can use Noether's theorem 
with 



L{p,p,t) = hp 2 -uj 2 p 2 -\) 
2 p z 



to prove that 



1 ~2 „2 

$ = _ 

-p z x 



- 2 + C 2 P- 2 + (px - pxf\ (8) 



2 

is a conserved quantity for Pinney's equation leading to 

x + uj 2 x = C 2 /x 3 . (9) 

The quantity C 2 is an arbitrary constant; choosing C 2 = 0, we reduce Pin- 
ney's solution to the parametric linear solution, while (||) turns into the 
Ermakov-Lewis invariant. 

If we write the invariant for two different solutions of the linear parametric 
equation, x\ and x 2 , while keeping the same p, and eliminate p in the resulting 
equations we get 

1 



p = ^\Jhx 2 2 + hx{ + 2x 1 x 2 [hl 2 - W 2 ] 1 / 2 , (10) 

where W — X\X2 — ^£-1^25 

and 1 1 and I 2 are constants. Thus, a general 
solution of Pinney's equation can be obtained if two solutions of the linear 
parametric equation can be found. (Since the Wronskian W is constant for 
two independent linear solutions, we can find that I\ — 1, J 2 = W 2 , and 



therefore (|T0|) turns into p = yxf + (1/W 2 )x 2 , which is the result given 
by Pinney in 1950). Moreover, one can see from (|7|) that two independent 
parametric solutions are x\ = pcos0, x 2 = psin0, where p is any solution of 



Pinney's equation. Then W — 1, and (pTOj) turns into 

p = py/h sin 2 + T 2 cos 2 + [hl 2 - lji/a sin 20 , = 1/p 2 . (11) 

This beautiful result obtained by Lutzky by means of Noether's theorem 
gives the general solution of Pinney's equation in terms of an arbitrary par- 
ticular solution of the same equation. Moreover, Lutzky suggested that this 
approach can be used to solve certain nonlinear dynamical systems once a 
conserved quantity containing an auxiliary function of a corresponding non- 
linear differential equation can be found. Even if the auxiliary equation is 
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nonlinear, sometimes it is simpler to solve than the original linear equation. 
In any case, one can establish useful relationships between the solutions of 
the two types of equations. 

In conclusion, we mention that Noether's method can be applied to the 
equation of parametric motion with first derivative ([I]); in this way one can re- 
produce the results of Eliezer and Gray of the previous chapter. The effective 
Lagrangian for (0) is given by L = |e F ^[x 2 — to 2 (t)x 2 ], where dF/dt = g(t). 

8. Possible generalizations of the Ermakov 
method. 

We have seen that there is a simple relationship between the solutions of the 
parametric oscillator 

X + LU 2 (t)x = , (1) 

and the solution of nonlinear differential equations of the Pinney type that 
differ from eq. (|1|) only in the nonlinear term. The equation of motion of 
a charged particle in some types of time-dependent magnetic fields can be 
written in the above form. Many time- dependent oscillating systems are 
governed by the same eq. (1). A conserved quantity for eq. ([I]) is 

lEL=\[{x 2 /p 2 ) + ( P X-pxf} , (2) 

where x(t) satisfies eq. ([J) and p(t) satisfies the auxiliary equation 

p + u 2 (t)p = l/p 3 . (3) 
Using eq. (P to eliminate ui 2 (t) in eq. (^|) we find 

p+(p/x)x = l/p 3 , (4) 

or 

xp — p'x = (d/ dt) (xp — px) = x/ p 3 . (5) 
Now, multiplying this equation by xp — px we can write 

(xp — px)(d/dt)(xp — px) = (xp — px)x/p 3 , (6) 
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or 



~(d/dt){xp - px) 2 = - l -{d/dt){x/p 2 ) , (7) 
and therefore we have the invariant 

lEL = \[{x 2 /p 2 ) + ( P X-pxf} , (8) 

where x is any solution of eq. (|1|) and p is any solution of eq. ([3]). 

A simple generalization of this result has been proposed by Ray and Reid 
in 1979 [|T2[]. Instead of (|3|) they considered the following equation 

p + u 2 (t)p=(l/xp 2 )f(x/p) , (9) 

where a; is a solution of eq. ([[]) and f(x/p) is an arbitrary function of x/p. 
If again we eliminate u 2 and we employ as a factor xp — px as a factor to 
obtain 

~{d/dt){xp - px) 2 = -(d/dt)(f)(x/p) , (10) 

where 

<t>(x/p) = 2 f P f(u)du . (11) 



From eq. fllCf ) we have the invariant 

I f = l -[ ( j ) {x/p) + {px-px) 2 } , (12) 

where x is a solution of eq. (P and p is a solution of eq. (§). For f = x/p we 
reobtain the invariant I el- The result ( |T2"D provides a connection between the 
solutions of the linear equation ([!]) with the solutions of an infinite number 
of nonlinear equations (|S[) by means of the invariant //. 

As an additional generalization, one can consider the following two equa- 
tions 

x + cu 2 (t)x = (l/px 2 )g(p/x) , (13) 

p + u 2 {t)p={l/xp 2 )f{x/p) , (14) 

where g and / are arbitrary functions of their arguments. Applying the same 
procedure to these equations one can find the invariant 

If, g = \l<P(x/p) + d{p/x) + {xp - px) 2 } , (15) 
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where 

<P(x/p) = 2f /P f(u)du, (16) 

6(p/x) = 2 g(u)du . (17) 

The expression ([T5]) is an invariant whenever x is a solution of eq. (|T^) and 



p is a solution of eq. flT4|). One should notice that the functions / and g 
are arbitrary, and therefore the invariant If t9 gives the connection between 
the solutions of many different differential equations. We can see that the 
Ermakov- Lewis invariant is merely a particular case of J/ >g with g — 0, f — 
x/p. 

In the cases g = 0, / = 0; g = 0, / = x/p; g = p/x, f = 0; and 
f — x/p, g — p/x the equations (fT3| ) and (|H] ) respectively are not coupled. 
In general, if we have found a solution for x, then the invariant J/ iS provides 
some information about the solution p. 

On the other hand, it is not known if the simple mechanical interpretation 
of Eliezer and Gray is also available for different choices of / and g. The 
simple proof of the existence of J/ j5 clarifies how such invariants can occur 
from pairs of differential equations. 



9. Geometrical angles and phases in the Er- 
makov problem. 

The quantum mechanical holonomic effect known as Berry's phase (BP) 
(1984) has been of much interest in the last fifteen years. In the simplest 
cases, it shows up when the time-dependent parameters of a system change 
adiabatically in time in the course of a closed trajectory in the parameter 
space. The wave function of the system gets, in addition to the common 
dynamical phase exp(— %h Jq E n (t)dt), a geometrical phase factor given by 

7n(c) =i [ dt(* n (X(t))\j t \* n (X(t))) , (1) 

because the parameters are slowly changing along the closed path c of the 
spatial parameter X(t) during the period T. \^/ n (X(t))) are the eigenfunc- 
tions of the instantaneous Hamiltonian H(X(t)). BP has a classical analogue 
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as an angular shift accumulated by the system when its dynamical variables 
are expressed in angle-action variables. This angular shift is known in the 
literature as Hannay's angle (Hannay 1985, Berry 1985). Various model sys- 
tems have been employed to calculate the BP and its classical analog. One 
of these systems is the generalized harmonic oscillator whoose H is given by 
(Berry 1985, Hannay 1985) 

H XYZ (p,q,t) = ^{X(t)q 2 + 2Y(t)qp + Z(t)p 2 } , (2) 

where the slow time- varying parameters are X(t), Y(t) y Z(t). 

Since Hxyz can be transformed into the if of a parametric oscillator, it 
follows that there should be a connection between the BP of the system with 
Hxyz and the Lewis phase for the parametric oscillator J/J. This problem 



has been first studied by Morales |L3| . Interestingly, the results appear to be 
exact although the system does not evolve adiabatically in time and goes to 
Berry's result in the adiabatic limit. 

Lewis and Riesenfeld [|7| have shown that for a quantum nonstationary 
system which is characterized by a Hamiltonian H(t) and a Hermitian in- 
variant I(t), the general solution of the Schroedinger equation 

ih^^- ^ H(tMq,t) , (3) 

is given by 

^ (q, t) = C n exp (ia n (t))V n (q,t) . (4) 

n 

ty n (q,t) are the eigenfunctions of the invariant 

i^ n (q,t) = \ n ^ n (q,t) , (5) 

where the eigenvalues are time-dependent, the coefficients C n are constants 
and the phases a n (t) are obtained from the equation 

hda n {t)/dt = (V n \ihd/dt - H(t)\ty n ) . (6) 

Using this result, Lewis and Riesenfeld obtained solutions for a quantum 
harmonic oscillator of parametric frequency characterized by the classical 
Hamiltonian 

H(t) = Ip 2 + ^ 2 (t)g 2 (7) 
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and the classical equation of motion 

q + n 2 (t)q 2 = 0, (8) 

where the points denote differentiation with respect to time. The matrix 
elements that are required to evaluate the BP are given by 

(* n \d/dt\* n ) = ±i(pp-?){n + ±) . (9) 

(*»|#(t)|*»> = \{(? + tf{t)p 2 + l/p 2 )(n + ~) , (10) 
where p(t) is a real number, satisfying the equation 

p + n 2 (t) P = i/p 3 . (ii) 

Substituting (Q) and ([U]) in (||) and integrating one gets 

a n (t)=-(n + h fdt/p 2 (t) . (12) 
2 Jo 

One can show this either by using (§) or ( |12|) and one can get the BP and 
Hannay's angle for the system of Hamiltonian Hxyz- For this system, the 
frequency that can be obtained from the Hamiltonian expressed in the action- 
angle variables, is given by 



to 



dH(I, X(t), Y(t), Z{t))/dl = (XZ - Y 2 ) 1/2 . (13) 



From (0) one can get the equations of motion for q and p and by eliminating 
p one can get the Newtonian equation of motion for q as follows 

q - (Z/Z)q + [XZ -Y 2 + (ZY - YZ)/Z]q = . (14) 

The term in q can be eliminated by introducing a new coordinate Q(t) given 
by (Berry 1985) 

q(t) = \Z(t)f 2 Q(t) . (15) 



Substituting fllBD in fljJD one gets 
Q+[XZ-Y 2 +(ZY-YZ)/Z+[1/2(Z/Z-Z 2 /Z 2 )-1/A(Z/Z) 2 ]]Q = , (16) 
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which corresponds to the equation of motion of an oscillator with parametri- 
cally forced frequency. Berry found Hannay's angle A8 by the WKB method 
in quantum mechanics, but it can also be obtained by means of (Q) or (|T^). 
Comparing @ with (|16|) we see that we can define Q 2 (t) as 

Q 2 {t) = XZ-Y 2 + {ZY-YZ)/Z+[1/A(Z/Z-Z 2 /Z 2 )-1/A(Z/Z) 2 } . (17) 

With this connection and employing (Jl|) and (|9|) we get 

7n(C) = ~(n+±) J* \pp - p 2 )dt , (18) 



where p(t) is the solution of (|Tl|) with Q 2 (t) given by (|17D. It is important 
to notice that ( |T8"D is exact even when the system does not evolve slowly in 
time. 

To compare with Berry's result one should take the adiabatic limit. For 
that we define an adiabaticity parameter e and a 'slow time' variable r 

f = £( r ) r = et , (19) 

in terms of which fi 2 (r) turns into 

n 2 (r) = e- 2 [XZ-Y 2 +e(Z'Y-Y'Z)/Z+e 2 [l/2(Z' /Z)'-1/A(Z' /Z) 2 }} , (20) 

where the primes indicate differentiation with respect to r. It has been shown 
by Lewis that in the adiabatic limit eq. (|20|) can be solved as a power series 
in e with the leading term given by 

p = ^ 1 / 2 (r). (21) 



If we plug this expression for p and its time derivatives in (|18D we could obtain 
the BP in the adiabatic limit. However, it is easy to calculate the Lewis 
phase first and then resting it from the dynamical phase — h~ l {^ n \H{t)\^! n ). 
Substituting (|2U|) and (F2T|) in (0) we get 

*<r) - ~(n + l) (I j\xZ - Y^W + \ /; Jglff^ r + 0(e)) 

(22) 
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The first term in the right hand side is the dynamical phase and the second 
and hiher order terms are associated with Berry's phase. Therefore we can 
write the BP as follows 

1 1 r T ZY — YZ 
l«(C) = --(n + -)j o z(xz _ Y2)I/2 dt . (23) 

Hannay's angle is obtained by using the correspondence principle in the form 
(Berry 1985, Hannay 1985) 

A# = -d ln /dn , (24) 



as 



which is the same result as that obtained by Berry (1985). 

In this way, it has been proved that if a time-dependent quadratic H can 
be transformed to the parametric form given by (0), then the Lewis phase can 
be used to calculate the BP and the Hannay's angle. Although we presented 
the particular case discussed by Morales, it is known that Lewis' approach 
for time-dependent systems is general. As a matter of fact, one can find more 
general cases in the literature. 



10. Application to minisuperspace Hamilto- 
nian cosmology. 

The formalism of Ermakov invariants can be a useful alternative to study 
the evolutionary and chaoticity problems of "quantum" canonical universes 
since these invariants are closely related to the Hamiltonian formulation. 
Moreover, as we have seen in the previous chapter, Ermakov's method is 
intimately related to geometrical angles and phases [p^j . Therefore, it seems 



natural to speak of Hannay's angle as well as of various types of topological 
phases at the cosmological level. 

The Hamiltonian formulation of the general relativity has been developed 



in the classical works of Dirac and Arnowitt, Deser and Misner (ADM) 



lj. When it was applied to the Bianchi homogeneous cosmological models 



it led to the so-called Hamiltonian cosmology pll. Its quantum counterpart, 
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the canonical quantum cosmology [JT7L is based on the canonical quantization 
methods and/or path integral procedures. These cosmologies are often used 
in heuristic studies of the very early universe, close to the Planck scale epoch 
t P « 10" 43 s. 

The most general models for homogeneous cosmologies are the Bianchi 
ones. In particular, those of class A of diagonal metric are at the same time 
the simplest from the point of view of quantizing them. 

Briefly, we can say that in the ADM formalism the metric of these models 
is of the form 

ds 2 = -dt 2 + e 2a ( t )(e 2 ^) ij o; i a; j , (1) 

where a(t) is a scalar function and /3y(t) is a diagonal matrix of dimension 
3, /3y = diag(x + v3y, x — V^y, — 2x), to 1 are 1-forms characterizing each of 
the Bianchi models and fulfilling the algebra du 1 = |Cj k c<jj A u k , where C] k 
are structure constants. 

The ADM action has the form 

Iadm = / (P x dx + P y dy + P a da - NWjdt, (2) 

where the Ps are the canonical moments, N is the lapse function and 

H± = e- 3a (-P 2 + P 2 + Pj + e 4Q V(x, y)) . (3) 

e 4ct V(x, y) = U(q M ) is the potential of the cosmological model under consid- 
eration. The Wheeler-DeWitt (WDW) equation can be obtained by canon- 
ical quantization, i.e., substituting P qM by P q M = — i<9 q M in eq. (3), where 
q M = (a, x, y). The factor ordering of e~ 3a with the operator P a is not unique. 
Hartle and Hawking [|18| suggested an almost general ordering of the follow- 
ing type 

_ e -(3-Q)a d a e- Qa d a = -e~ 3a d 2 a + Q e~ 3a d a , (4) 
where Q is any real constant. If Q = the WDW equation is 

□ tf-U(q")tf = 0, (5) 

Using the ansatz ^/(q^) = Ae ±$ one gets 

±An<S> + A[(V<S>) 2 -U} =0, (6) 

where □ = G""^, (V) 2 = + (£) 2 + y = diag(-l, 1, 1). 
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Employing the change of variable (a,x, y) — > (fli, 02, Pz), where 

Pi = a + x + v%, fa = a + x — V3y, (3% = a — 2x, (7) 

the ID character of some of the Bianchi models can be studied in a more 
direct way. 

Empty FRW (EFRW) universes for Q = 0. 

We now apply the Ermakov method to the simplest cosmological oscilla- 
tors which are the empty quantum universes of Friedmann- Robertson- Walker 
(EFRW) type. The results included in this chapter have been published re- 
cently When the Hartle-Hawking parameter is equal to zero (Q = 0), 
the WDW equation for the EFRW universe is 

fg -«-»»(«) =0, (8, 

where Q is the Misner time which is related to the volume of the universe 
V at a given cosmological epoch as fl = — IniV 1 ^ 3 ) [p0[| , k is the curvature 
parameter of the universe (1,0,-1 for closed, plane, open universes, respec- 
tively) and ^ is the wave function of the universe. The general solution is 
obtained as a linear superposition of modified Bessel functions of zero order 
in the case for which k = 1, = Ciio(f e" 2n ) + C 2 K {\e~ m ). If k = -1 

the solution will be a superposition of ordinary Bessel functions of zero order 
^iP) = CiJ (^e^ 2n ) + C 2 Y (^e~ 2Q ). C\ and C 2 are arbitrary superposition 
constants that we shall take for simplicity reasons equal Ci — C2 — 1. 

Eq. (H) can be transformed into the canonical equations of motion for a 
classical point particle of mass M = 1, generalized coordinate q = and 
moment p = ^ , and by considering Misner's time as a Hamiltonian time for 
which we shall keep the same notation. Thus, we can write 

dn = p (9) 
In = Ke q - (10) 

These equations describe the canonical motion of an inverted oscillator [k = 
1) and of a normal one (/t = —1), respectively [BTI], of Hamiltonian 

H<n) = £-Ke-"£. (11) 
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For the EFRW Hamiltonian the phase space functions T\ = T 2 = pq, y 

2 

T 3 = \ form a dynamical Lie algebra, i.e., 



H = J2h n (n)T n (p,q) , 



(12) 



which is closed with respect to the Poisson brackets {Ti, T 2 } = — 2Ti, {T 2 , T 3 } = 
-2T 3 , {Ti,T 3 } = -T 2 . The Hamiltonian EFRW Hamiltonian can be written 
now as 

if = 71 - /te 4n T 3 . 
The Ermakov invariant f is a function in the dynamical algebra 



I = ]Te r (ft)T r , 



(13) 
(14) 



and through the time invariance condition 

one is led to the following equations for the unknown functions e r (fl) 



o, 



(15) 



(16) 



where C^ m are the structure constants of the Lie algebra given above. Thus, 
we obtain 



ei = -2e 2 

-AO. 

e 2 = -«e e 1 - e 3 
e 3 = ~2Ke~ m t 2 . 



(17) 



The solution of this system of equations can be easily obtained by choosing 
ei = p 2 , which gives e 2 = — pp y e 3 = p 2 + -\, where p is the solution of 
Pinney's equation 

(18) 



ne 



P 3 



In terms of p(S~2) and using (6), the Ermakov invariant can be written as 
follows 



^kin Ipot 



(pp — pq) 2 q 2 p 4 r d ( ^ 



+ 



2p 2 



dVL \p 



+ 



1 /£ 

2 L, 



(19) 
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We have followed the calculations of Pinney and of Eliezer and Gray for p(Q) 
in terms of linear combinations of the aforementioned Bessel functions that 
satisfy the initial conditions as given by these authors. We have worked with 
the values A = 1, B = —1/W 2 y C = for Pinney's constants, where W 
is the Wronskian of the Bessel functions. We have also chosen an auxiliary 
angular moment of unit value (h = 1). Since I = h 2 /2,we have to obtain a 
constant value of one-half for the Ermakov invariant. We have checked this 
by plotting 1(0) for k = ±1 in fig. 1. 

Now we pass to the calculation of the angular variables. We first calcu- 
late the time-dependent generating function that allows us to go to the new 
canonical variables for which I is chosen as the new "moment" [|| 




(20) 



leading to 



S(q,I,e(Q)) = ^ + Jarcsin 



+ 



qy / 2Ip' 2 — q' 2 



(21) 



where we have put to zero the constant of integration. Then 




(22) 



Now, the canonical variables are 




(23) 



The dynamical angle is 




(24) 



while the geometrical angle (generalized Hannay angle) is 




(25) 
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The sum of A6 d and AO 9 is the total change of angle (Lewis' angle): 

AO* = 1^ \dtt' . (26) 
Jn Q p 2 

Plots of the angular quantities (24-26) for k — 1 are displayed in figs. 2,3, 
and 4, respectively. For k — — 1 we've got similar plots. 



i 1 1 1 1 i 1 j 1 1 t 1 1 1 1 t 



fig. 1: The Ermakov-Lewis invariant for Q = 0, k = ±1, h = 1. 
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fig. 2: The dynamical angle as a function of il. 
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fig. 4: The total angle as a function of fl. 



EFRW universes for Q ^ 0. 

We now apply the Ermakov procedure to the EFRW oscillators when Q ^ 0. 
These results have been reported at the Third Workshop of the Mexican 
Division of Gravitation and Mathematical Physics of the Mexican Physical 
Society |2|. 



It can be shown that the WDW equation for EFRW universes with Q 
taken as a free parameter is 

where, as previously, Q is Misner's time and k is the curvature index of the 
FRW universe; k — 1, 0, —1 for closed, flat, and open universes, respectively. 
For k = ±1 the general solution can be expressed in terms of Bessel functions 

= e~ 2an (c x I a (^e-™) + C 2 K a ^e- 2n )^j (28) 
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and 

*-(n) = (c^Z-e-™) + C 2 Y a { l -e~^ , (29) 

respectively, where a = Q/4. The case k — does not correspond to a 
parametric oscillator and will not be considered here. The eq. (29) can be 
turned into the canonical equations for a classical point particle of mass M = 
e Qn , generalized coordinate q — \& and moment p = e Qn ^J (i.e., of velocity 
v — 4f). Again, identifying Misner's time Q with the classical Hamiltonian 
time we obtain the equations of motion 

i=% = e"«S (30) 

= «<«-*»,. (3D 

as derived from the time-dependent Hamiltonian: 

H cl (n) = e' QnV - - /ce«- 4 > n ^ . (32) 

The Ermakov invariant /(fi) can be built algebraically to be a constant of 
motion. The result is 

m = (P 2 ) • y - (^P/i) • PQ + (e 2QQ P 2 + ^) • y ' (33) 

where p is the solution of Pinney's equation p + Qp — Ke~ 4n p = 6 ^f" . In 
terms of p± (f2) the Ermakov invariant for this class of EFRW universes reads 

Ww- 2 + 24- — l^*«-^«J 2^J ' 

(34) 

In the calculation of I^frw we nave use d linear combinations of Bessel func- 
tions that fulfill the initial conditions for p as explained in chapter 6 and 
which will be presented in some detail in a separate section of the present 
chapter. 

Calculating again the generating function S(q, I, Q) of the canonical trans- 
formations leading to the new momentum /, we obtain 



2 - 

S(q,I,0) = e Qn ^-^ + Jarcsin 
2 p 



VW 2 



+ 27 ' (35) 
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where the integration constant is again chosen to be zero. The new canonical 
variables are q\ = pV2Isin6 y pi = ^y-(^cos9 + e Qn pp sin 9*j. The angular 



quantities are: A6 d = fg(^)dQ' = - ^(e« n 'pp) + e^'p 2 ]dQ', 

= \ In Q [-^r(e Qn pp) — 2e Qn p 2 }dVt , for the dynamical and geometrical 

angles, respectively. Thus, the total angle will be A9 = Jq q 6 pi dVt . On the 
Misner time axis, going to — oo means going to the origin of the universe, 
whereas = means the present era. With these temporal limits for the 
cosmological evolution, one finds that the variation of the total angle A9 is 
basically the same as the Laplace transformation of 1/p 2 : A9 = —Li/ p 2(Q). 

The plots of the invariant and of the variations of the angular quantities 
are shown next, both for the closed EFRW universes as for the open ones. 




fig. 5: Jg FRW (f2) for Q = 3 and an initial singularity of unit auxiliary angular momentum. 



35 




fig. 6: The dynamical angle as a function of fl for a closed EFRW universe and Q = 1. 




fig. 7: The geometrical angle for the same case. 
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fig. 8: The total angle as a function of f2 for the same case. 




i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 t 



fig. 9: ^efrw(^) w ith Q = 1 for an initial singularity of auxiliary angular momentum excitation h = 2. 
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fig. 10: The dynamical angle as a function of for an open EFRW universe of Q = 1. 




fig. 11: The geometrical angle as a function of for the same case. 
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fig. 12: The total angle as a function of f2 for the same open case. 
Somewhat more complicated cosmological models 

We sketch now the Taub pure gravity model whose WDW equation reads 

d 2 m 



dn 2 



d 2 ^ 

+ O_ + e -V(0, = o, 



(1) 



where V(($) = |(e~ 8/3 — 4e _2/3 ). This equation can be separated in the vari- 
ables x\ = —4Q — 8/3 and x 2 = — 4S1 — 2/3. Thus one is led to the following 
pair of ID differential equations for which the Ermakov procedure is similar 
to the EFRW case 

' , .2 



d 2 t> 



i 1 Q d^Ti 



and 



dx 2 



d 2 ^T2 



12 dx! 

Qd^T2 



1 

144' 



n 







+ [u 



1 

—e 
9 



X2 



vJ/ T2 = . 



(2) 



(3) 



dx 2 3 dx 2 

where u;/2 is a separation constant. The solutions are \& T1 = ^Tai - 
e (-Q/u)xiZ iai (ie Xl / 2 / 6) and ^ T2 = ^ Ta2 = e^ Xa Z iaa (i2e B ^ 2 /3), respec 
tively, where «i = X 2 - (Q/12) 2 and a 2 = </4j 2 - (Q/3) 2 . 
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A more realistic case is that in which a scalar field of minimal coupling 
to the FRW minisuperspace metric is included. The WDW equation is 

[d 2 + Qd n Ke~^ + m 2 e- 6 VMft, <t>) = , (4) 

and can be written as a Schroedinger equation for a two-component wave 



function (see p3| ). This allows to think of squuezed cosmological states in 
the Ermakov framework p4| . For this, we shall use the following factorization 
of the invariant J = h(btf + ~), where b = {2h)~ 1 / 2 [ g + i(pp — e QcU pq)] 

and = (2/i) -1 / 2 [^ — i(pp — e Qcn pq)]. Q c is a fixed ordering parameter. 
Consider now a Misner reference oscillator of frequency u) Q corresponding 
to a given cosmological epoch for which one can introduce the standard 
factorization operators a = (2hu!o)~ 1 ^ 2 [u!oq + ip], a) = (2hu>o)~ 1 ' 2 [uJoQ — ip\. 
The connection between the two pairs a and b is b(Q) = p(Q)a + u(Q)a^ y 
&t(ft) = n*(n)rf + u*{Sl)a) , where //(Q) = (4a; )- 1/2 [p _1 - ie Q ^p + cu p] 
and z/(fi) = (4cj ) -1 / 2 [p -1 — ie QcQ p — u p\ satisfy the relationship \p{fi)\ 2 — 
|^(f2)| 2 = 1. The uncertainties can be calculated (Ag) 2 = ^\p— v\ 2 , (Ap) 2 = 
^l/i + z/| 2 , and (Ag)(Ap) = ||// + v\\fjL — v\ showing that in general these 
Ermakov states are not of minimum uncertainty [p4]| . 

The way one should do the linear combinations for the solutions 
of the linear differential equations. 



As it has been shown, in order to solve Pinney's equation one should first 
find the solutions to the equations of motion. Since these equations are lin- 
ear, we have chosen those combinations which satisfy the initial conditions 
of motion. According to the interpretation of Eliezer and Gray, the solu- 
tion of Pinney's equation is just the amplitude of the 2D auxilliary motion. 
Therefore, two of the three quadratic terms of the solution can be seen as 
the amplitudes along each of the axes, respectively. The third one is a mixed 
term (one can also eliminate it by diagonalizing the quadratic form in the 
square root). 

Let q(0) = a and g(0) = b be the initial conditions for the equation of 
motion. The solution can be written as x (t) = ax\ (t) + bx2 (£), and therefore 
the functions x\ y x 2 must satisfy the conditions X\ (0) = 1, X\ (0) = 0, 
x 2 (0) = 0, x 2 (0) = 1. If we take ipi an d ip 2 as a pair of linear independent 
solutions of the parametric equation, then we can build the functions xj 
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as linear combinations of ipi\ x\ = a\ipi + b\ip2- It is clear that the linear 
superpositions that satisfy the initial conditions will be: 



Xi 



W(0) 



4(0)^(^-^(0)^(0 



(i) 



x 2 



W(0) 



(2) 



where W(0) is the Wronskian of the functions ipi and ip2 evaluated at zero 
time parameter. The functions x\ are the correct ones that should enter the 
solution of Pinney's equation written in the form given by Eliezer and Gray. 
In this way, we have in the case of the cosmological models that have been 
discussed: 

(2z)* 



Xi 



V>i (1/2)Kq(z)-^ 2 (l/2)/ fi (z) 



(3) 



where 



x 2 



(2zY 



Kq (1/2)1 q(z) -Iq(1/2)Kq(z) 



1 



2Q. 



^'(1/2) = - 



^'(1/2) 



Q 



7fl(l/2) + / fl '(l/2) 



%Ka(l/2) + Ka' (1/2) 

Z 4 4 



(4) 



for closed EFRW, and similarly for the open EFRW models. 

The superposition coefficients we worked with are of the form a + = 
N K (l/2)/D + (l/2), b + = Nj(l/2)/D + (l/2), c + = -K(l/2)/ D + (l/2), d + = 
1(1/2)/ D + (l/2), where N K (l/2) = Kq +1 (1/2) - QKq(1/2), Nj(1/2) = 
Iq +1 (1/2)+QKq(1/2), and J D+(l/2) = Iq +1 (1/2)Kq(1/2)+Kq +1 (1/2)Iq(1/2) 

= -N Y (l/2)/D_(l/2),b_ = Nj(l/2)/D_(l/2), 
- J(l/2)/D_(l/2), where N Y (l/2) = V fi+1 (l/2)- 
QYq(1/2),Nj(1/2) = Js +1 (l/2)+gJs(l/2),and J D_(l/2) = Jq +1 (1/2)Yq(1/2)- 
Yq 1,(1 / 2) J q (1/2) for the open EFRW case. 



for the closed EFRW case; a. 
c_=y(l/2)/£L(l/2),d_ = 



41 



11. Application to physical optics. 



In order to study the Ermakov procedure within physical optics, our starting 
point will be the ID Helmholtz equation in the form given by Goyal et al 
p5f| and Delgado et al |26 



+ A0(a#(x) = , (1) 

that IS, cLS cL Sturm-Liouville equation for the set of eigenvalues A G R defin- 
ing the Helmholtz spectrum within a closed given interval [a,b] on the real 
line, where the nontrivial function if) turns to zero at the end points (Dirich- 
let boundary conditions). Eq. (1) occurs, for example, in the case of the 
transversal electric modes (TE) propagating in waveguides that have a con- 
tinuously varying refractive index in the x direction but are independent of 
y and z. Similar problems in acoustics can be treated along the same lines. 
The transformation of eq. (1) into the canonical equations of motion of a 
classical pointparticle is performed as follows. Let if){x) by any real solution 
of eq. (1). Define x = t, ip = q, and if) = p; then, eq. (1) turns into 

d <i /„x 
Tt = p (2) 

| = -A0(t)g, (3) 

with the boundary conditions q(a) = q(b) = 0. The corresponding classical 
Hamiltonian 

H(t) = £ + ty{t)£. (4) 

is similar to the previous cosmological case of Q = 0, if one identifies A = — k 
and (f> = e~ AQ . The procedure to find the Ermakov invariant follows step 
by step the cosmological case. In the phase space algebra we can write the 
invariant as 

i = X> r (*)r r , (5) 

r 

and applying 
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we get the system of equations for the coefficients p: r (t) 



/'i 

p3 



-2/i 2 

A0(t)/ii - /i 3 



(7) 



The solutions can be written in the conventional form by choosing Hi = p 2 , 
that gives /i2 = —pp and /i3 = p 2 + where p is a solution of the Pinney's 
equation of the form: p+\<p(t)p = with the Ermakov invariant of the well- 



2 ^ 2f? ■ Next, we calculate the generating function 



known form / — ^ pp ~^ 

of the canonical transformation for which / is the new momentum 

ri 



Thus, 



S(q,I,m) 



q 2 p 

S(q, I, fl(t)) — — h Jarcsin 

I p 



dqp(q ,I,fl(t)) . 



V2Ip 2 - q 2 



+ 



qVzTp 2 



2p 2 



(8) 



(9) 



where we have put to zero the integration constant. In this way we get 

OS . , q 

—— = arcsml 

ol 







in(- 



The new canonical variables are q\ = p\/2Ism9 and pi 
ppsinO). The dynamical angle is given by 



A6 a 



to 



2 dt' \p) 



dt 



whereas the geometrical angle is 

1 r* 

2 Jt 



(PP) - P 2 



dt 



(10) 

^(cosfl + 
p V 

(11) 
(12) 



For periodic parameters /!(£), with all the components of the same period T, 
the geometric angle is known as the nonadiabatic Hannay angle |2TJ that can 
be written as a function of p: 



A9 9 H 



fjdp. 



(13) 
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Now, in order to proceed with the quantization of the Ermakov problem, 
we turn q and p into operators, q y p = —ih-jp, but keeping the auxiliary 
function p as a real number. The Ermakov invariant is now a Hermitian 
constant operator 

d i = d i + klH)=V (14) 
dt dt ih v 1 v ' 

and the time-dependent Schrodinger equation for the Helmholtz Hamiltonian 
is 

ih^\mt)) = \(f + mW)\mt)) . (is) 

The problem now is to find the eigenvalues of I 

I\i> n (q,t) = K n \if) n (q,t)) (16) 
and also to write the explicit form of the general solution of eq. (|15|) 

mt) = J2c n e tan{t) Mq,t) (17) 

n 

where C n are superposition constants, ip n are (orthonormalized) eigenfunc- 
tions of I, and the phases a n (t) are the Lewis phases [[13], Q that can be 
found from the equation 

S ^ = W 4"^ W - (18) 
The crucial point in the Ermakov quantum problem is to perform a unitary 
transformation in such a way as to get time-independent eigenvalues for the 
new Ermakov invariant I = UIW. It is easy to obtain the required unitary 

a - ■ -2 '"v 2 "2 

transformation: [/ = exp[— f^y]- The new invariant will be J = ^y- + 

^-if. The eigenfunctions are oc e~2B H n (9/y/h), where H n are the Hermite 
polynomials, 9 = -, and the eigenvalues are n n = h(n + |). Thus, one can 
write the eigenfunctions ip n as follows 

^ « ph exp (^^ 2 ) exp ( - 4^) Hn ^- P ) ■ (19) 

The factor 1 / p 1 / 2 has been introduced for normalization reasons. Using these 
functions and doing simple calculations one can find the geometrical phase 

1 , 1, 



at = — in + - 
2 V 2' 



to 



(pp) - p 



dt . (20) 
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The cyclic (nonadiabatic) Berry's phase [^7] is 



a 



B,n 



\)j c Pdp. 



(21) 



The results obviously depend on the explicit form of p which in turn depends 
on the explicit form of <fi. 

One can find that a good adiabatic parameter is the inverse of the square 
root of the Helmoltz eigenvalues, -7=, with a slow "time" variable r = -^=b. 
The adiabatic approximation has been studied in detail by Lewis 0. If the 
Helmholtz Hamiltonian is written down as 



# (t) = + <Kt)g 2 ] 



(22) 



then Pinney's equation is 



-p + (f)(t)p 



1 



while the Ermakov invariant becomes a 1/vA- dependent function 



(23) 



J(1/VA) 



(pp-pq/VX) 2 , q 



+ 



2p 2 



(24) 



In the adiabatic approximation, Lewis obtained the general Pinney solu- 
tion in terms of the linear independent solutions / and g of the equation of 
motion jq+fl 2 (t)q = for the classical oscillator (see eq. (45) in |J). Among 
the examples given by Lewis, it is Q(t) = bt m / 2 , m ^ —2, b = constant which 
is directly related to a realistic dielectric of a waveguide since it corresponds 
to a power-law index profile (n(x) oc x m ^ 2 ). For this case, Lewis obtained a 
simple formula for p of 0(1) order in 1/\^X 



where H 



(i) 



and Hf 



7i 



7 2 7T 



Vx 



(m + 2) 



t*[H$\y)Hf\y)] 



(25) 



2bV\ t f +1 d 
(m+2) 1 ' clIiU '! 



are Hankel functions of order (3 = l/(m + 2), y = 
±1, 72 = ±1. An even more useful technological 
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application might be the following proposal of Lewis: 



m 



■in 



±1.±2. 



leading to 



2n+l' 



Pn = 7i7fHt^|G(M/\/A)| : 



where 



G(t, l/VX) 



fc=0 



-r 



n + k)\ i 1/y/X 



■( 



k\(n - k)\ ^2ib(2n + 1] 



t (2n+l) 



(26) 



(27) 



One gets p as a polynomial in the square of the adiabatic parameter, i.e., 
A -1 , of infinite radius of convergence. The topological quantities (angles 
and phases) can be calculated by substituting the explicit form of Pinney 's 
function in the corresponding formulas. Lewis || found a recursive formula 
in 1/A of order 1/A 3 that can be used for any type of index profile. The 
recurrence relationship is 



P = Po + Pi/A + p 2 /A 2 + p 3 /A 3 + ... 



(28) 



where po = f2 1 / 2 



-1/4, 



x)\ for the other coefficients pi see the appendix 



in p. The main contribution to the topological quantities are given by po- 
In the case of a power-law index profile, the geometric angle is 



A9 9 



m 



Ab(m + 2) 



- (f+1) _,-(?+!) 



(29) 



and a similar formula can be written for the geometric quantum phase. For 
periodic indices, one can write the Hannay angle and Berry's phase accord- 
ing to their cyclic integral expressions. Finally, we notice that the choice 
4>{x) = $(x) + ^§^y, which corresponds to nonlinear waveguides, leads to 
more general time-dependent Hamiltonians that have been discussed in the 
Ermakov perspective by Maamache p8| . 

We have presented in a formal way the application of the Ermakov approach 
to ID Helmholtz problems. For more detailes one can look in a recent work 
by Rosu and Romero [^9] . 
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12. Conclusions. 



As one could see from the examples we discussed in this work, the Ermakov- 
Lewis quadratic invariants are an important method of research for para- 
metric oscillator problems. They are helpful for better understanding this 
widespreaded class of phenomena with applications in many areas of physics. 
One can also say that the Ermakov approach gives a connection between 
the linear physics of parametric oscillators and the corresponding nonlinear 
physics. 

The cosmological applications of the classical Ermakov procedure we pre- 
sented herein are based on a classical particle representation of the WDW 
equation for the EFRW models. We also notice that the Ermakov invariant is 
equivalent to the Courant- Snyder invariant of use in the accelerator physics 
[ p0|1 , allowing an analogy between the physics of beams and the cosmological 
evolution as suggested by Rosu and Socorro |H[ . 



We end up with a possible interpretation of the Ermakov invariant within 
the empty minisuperspace cosmology. If one performs an expansion of the 
invariant in a power series in the adiabatic parameter, the principal term 
which defines the adiabatic regime gives the number of adiabatic "quanta" 
and there were authors who gave classical descriptions of the cosmological 



particle production in such terms [32]. On the other hand, the Eliezer-Gray 



interpretation as an angular momentum of the 2D auxiliary motion allows one 
to say that for EFRW minisuperspace models, the Ermakov invariant gives 
the number of adiabatic excitations of the auxiliary angular momentum with 
which the universe is created at the initial singularity. 
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Appendix A: Calculation of the integral of /. 

The phase space integral of I in chapter 5 can be calculated from the formula 
(15a) in the paper of Lewis || 

1 r 27r dXi 
2n Jo dtp 

where X% and X 2 represent the functional dependences of q and p, respec- 
tively, in terms of the nice variables z\ and (p, which have been given by 
Lewis in the formulas (38) of the same paper as follows 

Xl = ± F 1 fi[l+tan 2 ( V 9-F 2 )] 1 /2 ^ 

and 

_ Zl [e*±* + j,tui(<p - F 2 )] 

2 F^fl+tan 2 ^-^)] 1 ^ ' [6) 

where F\ and F2 are two arbitrary functions of time. Thus, 

V = ± ^ [1 + tan2( ^ ~ F2 )] _3/2 (- 1 /2)2tan(^ - F 2 )sec% - F 2 ) . (4) 
We have the following integral 

Now, employing 

s = tan 2 (y? — F 2 ), ds = 2tan(<y? — F 2 )sec 2 (y? — F 2 )d(p (6) 

one gets 

r or f - 

;i + s ) 2 P 2 J + 



1 « 77^ + 3 / 77—^. • (7) 



Therefore 

dlnp 



z 2 

I 



2?rF 1 2 2 



— e- 



dt (1 + 



-— ^ + - (-s" 1 / 2 + tan-yi 
1 + s ) cr v 



(8) 
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Going back to the <p variable and taking into account the corresponding 
and 2n limits one gets 

(9) 



which is the result obtained by Lewis. The common form of I can be obtained 
by going back to the (q,p) variables. 



q = ' ' 

{p + i/p)a + + {p-i/p)a 



Appendix B: Calculation of the expectation 
value of H in eigenstates of /. 

From the formulas (12) in chapter 5 for the raising and lowering operators 
one gets 

(1) 
(2) 

(3) 



V 



V2 



Performing simple calculations, one gets 
H = f(p)a +2 + r(p)a 2 + 1 - 



* 2 i i 2 2 

p + — + UJ P 



(27) 



where f(p) = p 2 + 2ip/p — 1/p 2 + uj 2 p 2 . Thus, 



(n\H\n) = (n\- 
2 



p 2 + — + uj 2 P 2 

p 2 



I\n) 



(4) 



from which eq. (16) in chapter 5 is obvious. 
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